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Abstract 



Quantum mechanical interference of wave functions leads to some dif- 
ficulties if a probability density is considered as a source of gravity. We 
^ I show that an introduction of a quantum energy-momentum tensor as a 

source term in Einstein equations can be consistent with general relativity 
if the gravitational waves are quantized. 
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Q>^ ; 1 Introduction 



The conceptual difficulties with an introduction of quantum effects into Einstein 
gravity are well-known. It would be a pessimistic attitude to think that the 
inconsistencies can be avoided only in a complete quantum theory of gravity 
and matter . It has been shown that some semiclassical methods to combine 
classical gravity and quantum mechanics contradict experiments and intuition 
^ . • Page and Geilker |^ suggested that decoherence could help to solve these 

?—( ' difficulties. An interaction with the environment as a source of the decoherence 



has been discussed in l|]-|§l. In ref.g Einstein gravity has been discussed and 
a phenomenological stochastic coupling has been introduced in order to achieve 
the decoherence. 

Our aim in this paper is to discuss quantum mechanics combined with the 
linearized quantum gravity as a consistent approximation to a hypothetical com- 
plete quantum theory. First, let us outline a scenario implying a physical rel- 
evance of the quantized gravity. We consider a model of a gravitational field 
interacting with the quantum complex scalar field. If these fields result from a 
quantization of the classical field theory then they should satisfy the operator 
equations of motion [n = SttG, where G is the Newton constant) 

DgV = M'^c^iP (1) 






where g^i, is the Riemannian metric, R^i, is the Ricci tensor constructed from 
the ChristofFel symbols TP and the wave operator Ug is 

n5=3^"(5^+.9"%a0)a. (3) 

Subsequently we can consider a reduction of the quantum field theory of the 
scalar field to quantum mechanics . For this purpose we restrict the field equa- 
tions to a fixed iV-particle sector. Let '4>ci be a solution of the Klein-Gordon 
equation. By means of the asymptotic fields we can define the one particle state 
IV'c;) • Then, for the matrix elements {0\ip{x)\^ci) we should obtain from eq.(l) 
approximately the Klein-Gordon equation in the metric g . Taking the expecta- 
tion value of both sides of cq.(2) in the state |V'c/) we obtain the (quantum) Ein- 
stein equations determined by the energy-momentum tensor {ipci\Tp,u{''P)\''Pci) of 
the quantum field theory. 

At this stage it is also useful to see the non-relativistic limit of eqs.(l)-(2). 
The non-relativistic limit is simple if g^^ = and g^^ is a slowly varying function 
of time. Then, we write 

(7°" = -l + y/c2 (4) 

and 

tp = exp( — ^^ojV 

As a result we obtain an equation for ip ( where A_^ is the thrcedimcnsional 
Laplacian on a Riemannian manifold M) 

ihdti^^{-^A^^-V)i^ (5) 

i.e., the conventional Schrodinger equation in the gravitational field. Then, in 
the limit c — > oo we obtain an equation for V 

We consider a semiclassical approximation to the solution of the Klein- 
Gordon equation on a manifold. The semiclassical approximation is determined 
by the Hamiltonian 

n = ^5^"^('Z)PmP. (6) 

A solution of the Hamiltonian equations of motion determines a solution of the 
Klein-Gordon equation in the leading order of h . The Hamilton- Jacobi function 

W{x) = f p^.dq'^ 



is computed on the trajectory starting from x and ending in y. We have in the 
semiclassical approximation 

V'(x) =exp(jy%^dg^)(det(^))^0(q.(a;)) (7) 

where qr is the trajectory starting from a; at t = and ending in y at r. We 
have to exchange r in favor of xq in eq.(7) . 

The Hamiltonian equations have a constant of motion which we express in 
the form 

For an approximate solution of the Hamihonian equations it wiU be useful to 
introduce the frames e 

e'aV'^'el = r (8) 

where 77 — (1, 1, 1, —1) is the Minkowski metric. We express e by a 

e^a = 77^c(expa)^ 

Then, in the lowest order 

e^a = ?7Ma + nv^cal (9) 

Geodesic equations (which are equivalent to Hamiltonian equations) can be 
expressed in the form 

d{e^a^) - (10) 

de^a+T''^ae.adq"^0 (11) 

We solve these equations by iteration . So, we have from eq.(lO) 

-^ = -e^'^PalM (12) 

(XT 

In the zeroth order in a the trajectory is a straight line q = x ~ Pt/M . We 
obtain a relation between t and xq which can be applied in order to eliminate 
T from W in eq.(7) 

xo-yo^ tPq/M 

2 Superposition principle and gravitational fluc- 
tuations 

We suggest that a quantization of gravitational waves may help to reconcile 
classical Einstein gravity with quantum mechanics. The common way to intro- 
duce a quantum matter into the classical gravity B is a replacement of T^^, by 
{'^\Tfj,^\ip) ■ However, we encounter the Schrodinger cat problem with such an 



extension of classical gravity. We obtain a superposition of states with a certain 
density determined by jV'P- In the classical gravity the metric is determined by 
the energy distribution (possibly random) , but not by the probability density 
of an interference of states of different energies [0 . 

We consider a perturbation expansion of eqs.(l)-(2) g — g'^'^'' + g^"?^ with 
the approximation of a weak classical field g'^K Then, □g^'^-' = 0. Hence, the 
quantum gravitons are approximately moving as waves in a flat background. 
Then, expanding in g^'^'"' 

n^<,,9(f - 2i?(«)^5^:') + 2Ri%J^''^^^ = ^\Mt) + Mtr (13) 

where we take a superposition ipi + ip2 of solutions of eq.(l). The observed 
value of the gravitational field is an average over the quantum fluctuations of 
the metric 

(n^<,,5(f - 2i?('')^5^:') + 2R^^l,g(^^^n = ^{\Mt) + Mt)\') (14) 

In the linear approximation the l.h.s. of eq.(14) is ng^'^^\ In the semiclas- 
sical approximation for the solution of the Schrodinger equation we obtain a 
superposition of solutions of eq.(7) 

A +^2^ exp( J W^(i) (t))0(g(i) it)) + exp( Jm^(2) (t))0(g(2) (i)) 

where g^'^^ (i) are the classical trajectories. At the point x the probability density 
(!?/'( (x)p) is equal to the diagonal part of the density matrix 

Pt^{\A){A\) (15) 

We assume that the quantum gravity consists of asymptotic flelds a'^ de- 
scribing gravitational radiation and eventually some other degrees of freedom 
which do not have asymptotic fields. We quantize at the moment only the grav- 
itational waves. In the linear approximation the metric as well as the tetrad e 
(a in the linear approximation) satisfy the wave equation 

We quantize the gravitational waves expanding a in the momentum space 

a^(a;) = V4^\/Wc^(27i-)~5 / dk|k|-^ exp(ikx) 



Ec=i,2 {^aiC, k)C{C, k) exp(-z|k|a;o) + S'aiC, k)C{C, k)+ exp(z|k|xo)) 
where 

[C(c,k),c(c',k')+] = 5cc"5(k-k') 

are the creation and annihilation operators. 



(16) 



The Hamiltonian is 



Hr = I dkJ2c\k\c{c,k)+c{c,k) 



c 

If the gravitational radiation is in equihbrium with hght and matter then it 
should be described by the Gibbs distribution 

P0 = Z-^exp{-f3HR) 

where 4 = KT , K is the Boltzmann constant and T denotes the temperature. 
In cosmological models the Gibbs distribution is believed to be correct when 
applied to the primordial gravitons (present at the earliest stages of the big bang) 
which had time to reach an equilibrium with other particles [Q . We suggest that 
these relict gravitons now reach the Earth and impose the classical behavior of 
large quantum systems. During an expansion of the Universe some gravitons are 
continuously being created as a result of a time-dependent gravitational field. 
These gravitons will not reach any equilibrium with the primordial ones. They 
have another energy distribution. The probability of graviton production is 
large at small wave number k. Hence, the Planck distribution may be modified 
for small wave numbers k [pj. For these reasons we consider a more general 
density matrix p(Hr) as a function of the graviton energy Hr. We introduce a 
parameter 1/6 (with the dimension of the inverse of the energy ) as an energy 
cutoff . We could represent p by a Fourier-Laplace transform of the Gibbs 
distribution p^ (fi may be complex) 

P= dp-f{p)p^ (17) 

Then, for a computation of expectation values we can use the methods applied 
for the Gibbs state , e.g., by a direct computation through an expansion in the 
number states we find in the Gibbs state 

/^^^(?ic|k|) ^ {C+ik)Cik))p = (exp(/3c?i|k|) - l)"' (18) 

We can see that effectively 1//3 plays the role of the energy cutoff in the Planck 
distribution because f?'"{hkc) decays exponentially fast when kch > 1//3. We 
shall often identify b with (3 in our discussion. 

The correlation functions can be computed using eqs.(17) and (18) 



Tr(<(i,x)<,(0,x')/5 1 =G^(x,x';i)^-, = 
2^ I d^^cc' ^co5(k(x - x')) ^ (19) 

(i -I- ft{hc\k\))cos{c\k\t) ~ ism(c|k|t) 



where we denoted (a representation of the tensor S and its properties are dis- 
cussed in Q) 

c 

In eq.(19) fb is the graviton distribution. Our results for small time and space 
separations do not depend essentially on the form of fi, if 

Mk) = f{bk) 

and if / decays sufficiently fast, e.g., |/(fc)| < Ak~^ for a large k. For a large 
time and large space separations the results depend on the singularity of f{k) 
at fc = (for the Planck distribution / « k~^). The distributions derived in 
inflationary models [n2||jl3 behave powerlike in some intervals ,e.g. , 

/(fc) = 

if fc > 1 and 

/(fc) = k-'^ (20) 

if < fc < 1 . Such a distribution leads to similar conclusions as the Planck 
distribution. However, an introduction of the infrared cutoff /(fc) = if fc > e 
would destroy the decoherencc at a sufRciently large time. 

An expectation value in the ground state x of the free gravitational field is a 
special case of eq.(19) corresponding to the limit /3 — > oo (under the assumption 

/oo(fc) = 0) 

Goo(x,i;x',0) =< x|a(t,x)a(0,x')|x >= 



KC 
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J dk-rh cos(k(x — x')) exp(— ic|k|t) 



(21) 



■|k| 

Comparing eqs.(19) and (21) it can be seen that the first term on the r.h.s. of 
eq.(19) describes the zero point density (vacuum fluctuations) whereas the sec- 
ond one comes from the thermal gravitons in equilibrium with the environment. 
In general, the vacuum fluctuations cannot be neglected. After a renormaliza- 
tion they contribute to measurable effects. However, it can be shown |jll| that 
renormalized vacuum fluctuations give a negligible contribution to the decoher- 
ence. They are small in comparison to the black body radiation at moderate 
temperatures. Moreover, the renormalized vacuum fluctuation part decreases 
to zero when the time becomes large. We subtract the vacuum fluctuations in 
eq.(19). After this subtraction the correlation function becomes real. We can 
associate a real random field a with such a correlation function (the field will 
be Gaussian in a linear approximation to gravity) 

where 

GZix - x') = ?iKC-3(27r)-2 / dk\k\-^S^Ak) cos(k(x - x')) cos((xo - a;[,)|k|)/f,(c|k|;i) 

(22) 



In spherical coordinates dk — 2Trddsm9dkk^ . Hence, Gth can be expressed in 
the form ( we skip the tensor S) 

Gth = 2hKC^^\x — x'l^^TT^^ Jg dksin{k\x — x'|) cos{ckt) fb{fick) 

3 Decoherent effect of gravitons 

It is natural to associate gravitons with the decoherence. Gravitons interact 
with all particles. Hence, their decoherence effect would be universal. We 
consider the Einstein gravity for weak fields . We define the partial density 
matrix (averaged over the gravitons) 

pt(x,x') =Tr^((x|/3(ifH)exp(-i^//fl)|V't)(Vt|exp(z|iJH)|x')) (23) 

The trace in eq.(23) can be obtained as an expectation value over the Gaussian 
random field a (or calculated in the operator formalism by means of the time- 
ordered products in the Fock space) 

(expiaJ) — exp(— -JGJ) 

where the Green functions G depend on the state under consideration. In partic- 
ular, in the thermal state with subtracted vacuum fluctuations G — > Gth ■ If the 
vacuum fluctuations were to be taken into account then we would need to make 
the replacement Gth —^ G/s = Gth + Goo and subsequently Goo —^ Gp = i/^F 
(in the notation of Bjorken and Drell[l^). However, the part J Gpdydy con- 
tains infinities when the paths intersect . After a renormalization the remaining 
expression gives a negligible contribution to the decoherence |ll[ . 
We consider an initial state 

^(x) = exp(iP(i)x/a)0(x) + exp(iP(2)x/?i)0(x) (24) 

We solve the geodesic equation (12) with the initial position x and the initial 
momenta P^'^' by iteration. Then, till the first order in a 



W ^^- ±pik) _L f a{x~ —p('^))pi'^)dT 



^" M' M 



The Green's function on the trajectory (till the zeroth order in a) is 

T — S T — S 

{a {y (r, x)) a {y (s, x'))) = G(x - x' ]f^P' a^o - 4 J^^o) 

In our approximation the time evolution is 

V'(x) ^ V*(x) = exp(*p(i)2/fV?i)'/'(yrVexp(zP(2)yfV?i)0(yf') 



where yl (x) = (yj {yi)^ct) and 



p(k) ^ /p(fe) /^j2^2 + (p(fe) 



Then, a calculation of the trace (23) leads to the formula ( we assume that (j) is 
a slowly varying function, hence in its argument it is sufficient to calculate y^'^-' 
in the lowest order in a) 

(I Vt(x) P) = \4>{^ J,P(1))P + |</.(X - ^P(2))12 , 



(x - i^P(2))0(x - ^P(i)) exp (A(P(2) _ P(i)) 



+0(x - i^P(i))0(x - j^P(2)) exp ( - i(p(2) - P(i))), 
exp(^/>Wp(i)G*.(^pW - ^P(2),p(^)./Af -Pf r/M)P(2)p(2)d.dr 
-2w/oP^'^P*'^G,.(^pW-^pW,p(^)(.-r)/M)pWpWd.dr 



2M2n^ 


Jo ^ 






Pi^' + 


pi^' 


+ 


(12' 
Pt 



(25) 
The notation PP is only symbolic, it means that we must sum the indices of P 
with the indices of the Green's function Gth ■ 

If the off-diagonal terms vanish then the density of a superposition of wave 
functions is a sum of the densities. Classically the packets move apart and 
each of them is independently a source of a gravitational field. Without the 
decoherence there would be a superposition of the states of the packets. 

It is easy to estimate the behavior of eq.(25) for a small time and either 
p(i) II p(2) or p(i) _L p(2). Let us denote in eq.(25) 

|p(i2)| ^ 2exp(-5i2(P))|</'(x- ^P(i))0(x- 1:P(2))| 

where 

S{P) = P'^Pf^Sl'i^PcP^ 

Then, for a small time we can set s = r = 0. There remains (under the 
assumption fb{u) = f{bu) ) 

5i2(P) = t2M-2n-2|(P(l))2 _ (P(2))2|2G,^(0,0) 

= Ai2M-2a-3|(p(l))2 _ (p(2))2|2^^-5^-2 J- duufi,{u/b) (26) 

= i|(P(l))2 - (p(2))2|(^)2y |(p(l))2 _ (p(2))2|;,-2^2^^ 

In eq.(26) A and A are constants of order 1, IdB = ^cb = Icbmc^ will be called 
de Broglie length because if f3 = b then I^b = ^cf3 is the wave length of a particle 
in a medium at temperature T, Ic = — is the Compton length and m is the 
electron mass as our basic mass unit. Then, ^ is particle's wave length at the 
momentum P , LpL = ^/finjc? is the Planck length. 



We can see that the decoherence is determined by the relation of the parti- 
cle's wave length and the length of particle's trajectory to the Planck's length 
and de Broglie's length. For a large time the calculations are more involved. 
Let us denote by 5*12 the expression Si2{P) without the four linear momenta. 
Then, we obtain a rather complicated formula for 6*12 (we omit the tensors S) 

( (P(2))' (kP(2)/Af + cky^ (1 - cos (ikP(2)/Af + tck)) 

+ (P(2))' (kP(2)/M - cky^ (1 - cos (tkP(2)/Af - tck)) 

+ (P(i))' (kP(i)/M + cky^ (1 - cos (tkP(i)/Af + tck)) 
+ (P(i))' (kP(i)/A/ - cky^ (1 - cos (tkP(i)/Af - tck)) 

+ 2(PWp(2))(p(l)k)(P(2)k)/M2(c2fc2_(p(l)k)VM2)"'(c2fc2-(p(2)k)'/M') 

(l-cos(ik(pW-P(2))/M)) 

-(P(i))2(p(2))2(c/c + p(i)k/M)-i(c/c + P(2)k/M)-i(l - cos(cA:t + tP(i)k/Af)) 

-(P(i))2(p(2))2(c/c _ p(i)k/Af)-i(c/c - P(2)k/Af)-i(l - cos{ckt - tP(i)k/A//)) 

-(P(l))2(P(2))2(c/c + P(l)k/Af)-I(c/C + P(2)k/TO)-l(l - COs{ckt + tP'^'^^k/M)) 

-(P(i))2(p(2))2(c/c - p(i)k/M)-i(c/c - P(2)k/Af)-i(l - cos(cA:t - tP(2)k/M))^ 

(27) ' 
It can be seen from eq.(27) that if /(O) = then owing to the Lebesgue theorem 
Si2{t) is bounded in t. If f{k) is singular at fc = then 5'i2(i) grows for a large 
t. In order to investigate this case in more detail let us consider first the Planck 
distribution fo^ . We neglect P/Mc in the denominator of eq.(27) and consider 
either P'^^ |1 p(2) or P^^^ _L p(2). Then, we obtain 

Si2iP) = M^^^li^^'^? iP^'^ri'Io" fieMPnck) - l)-i (1 - cos(teA:)) 
= jj^l^MP'^'^r-{P^'^??ctydaf^^dk{eMPnck)-iy^sin{ackt) 
= M4^|^l(P(^))^-(P«)^P/o'da(^coth(^)-5i,) 

= mW^IfKP^^^)^ - (P«)^Pln (^sinh(ll)) 

« _£^L2^^?i-2|(p(2))2 _ (p(l))2|2^-2^-2 

(28) 
for a large t. Here, the formula 3.911 of Gradshtein and Ryzhik M has been 
applied 

du sm(au) (^ exp(/3u) - 1 j = — coth( — ) - — 







For the inflationary distribution (20) f^ in eq.(28) is replaced by (j^)'^ . 

As a result of the decoherence the gravitational field is determined by the 
density distribution of two separate quantum particles. Hence, we have ob- 
tained a classical addition of probabilities instead of the quantum addition of 
amplitudes. This conclusion follows from the results (25)-(28) which should be 
inserted into eq.(14). 



Finally, let us discuss numerical estimates on some expressions in this paper. 
First, let us assume that the graviton distribution is determined by the Planck 
distribution (as suggested by Weinberg m). The decoherence time depends on 
the ratio of the Planck length to the de Broglie wave length (as could have been 
expected because we have two universal length scales at finite temperatures). 
To see the meaning of I^b let us write KT — mv^ then Ij^g — Ic^ (the 
Compton length multiplied by the ratio of the rest energy to the kinetic energy). 
For numerical estimates it is useful to introduce the Planck temperature Tp^ 
determined by Lpl — IdB, then /3^c^ = k/H and Lpl/IcIB = T/Tpl (note that 
TpL = 1.3 X 10"^^ Kelvin). Instead of the de Broglie length we could use the 
temperature independent length scale : the Compton length Ic (choosing the 
electron mass m as a mass unit). The Compton length comes in a natural 
way because the ratio of the gravitational force to the Coulomb force between 
electrons is 

Fgrav/Fcoul = -^(Lpl/Ic) ~ T^IO 



, .„ . - \2 _ ^ in-46 

grav / 
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In particular, we can write for a small time 

The exponential is of the order (10~^^ — )^(-)^. We can see that we need the 
number of particles TV = — « 10^'^ if the decoherence is to be visible . 

If instead of the Planck distribution we have the (simplified) inflationary 
distribution (20) then the behavior of 5i2 is determined by the integral (with a 
certain constant B) 

rb 



Si2 = B I dkk-"-^ (1 - cos {kct)) 
Jo 



Under the assumption < a < 2 this integral behaves as ( j^)*^ for a large t and 
as (j^)^ for a small t (we use the notation IdB — bch). In inflationary models a 
is different in various frequency ranges and for the lowest frequency range ct = 2. 
For such a a the integral 5*12 is infrared divergent. At the large wave length \/k 
there can be a sharp cutoff or a continuous change of the behavior of / leading 
to the finite integral for 6*12. This is a sensitive problem in the spectral theory of 
gravitational waves |14[ p5| | concerning waves with length larger than the actual 
Hubble horizon length. It is interesting that in this way the decoherence is 
associated with the prospective evolution of the Universe. 
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